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Abstract 

^ . Spin-orbit and spin-spin effects in the gravitational interaction are treated in a close 

analogy with the fine and hyperfine interactions in atoms. The proper definition of the 
cener-of-mass coordinate is discussed. The technique developed is applied then to the 
gravitational radiation of compact binary stars. Our result for the spin-orbit correction 
differs from that obtained by other authors. New effects possible for the motion of a 
spinning particle in a gravitational field are pointed out. The corresponding corrections, 
nonlinear in spin, are in principle of the same order of magnitude as the ordinary spin- 
spin interaction. 
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1. It is expected that in few years the gravitational radiation from coalescing binary stars 
will be observed by laser interferometer systems LIGO and VIRGO. Its successful detection 
depends crucially on the accurate theoretical prediction of the exact form of the signal. In 
this way the observed effect becomes sensitive to the relativistic corrections of the c -2 , c -3 
and c -4 orders to the motion of a binary system and to the radiation intensity. In particular, 
the spin-orbit interaction becomes essential, and for two extreme Kerr black holes even the 
spin-spin one [jj. 

Some years ago it was noticed that the general relativity can accomodate in a natural 
way a specific gravitational magnetic moment coupling || (see also 0). The starting point 
of the present work was the observation that the spin self-interaction arising in this way is 
of the same order of magnitude as the spin-spin interaction, and therefore in principle its 
existence can be checked in the gravitational- wave experiments. 

However, in the course of the investigation, when trying to rederive previous calculations 
related to the spin effects in the gravitational radiation of binary stars, we came to the results 
which differ from those of Refs. [jl], |]] as concerns the spin-orbit contributions. The origin 
of this discrepancy can be traced back to, what is to our belief, a long-standing confusion 
concerning the definition of the centre of mass in the case when spin is taken into account. 
The problem is quite instructive and amusing by itself, and on the other hand, the spin-orbit 
correction is the leading one among spin effects. That is why we would like to start our 
discussion with this subject. 

2. The spin-orbit and spin-spin interactions in the two-body problem can be immediately 
obtained in fact from the well-known results for the limiting case when one of the bodies (say, 
2) is very heavy (see, e.g., book ||). In this limit we have the usual spin-orbit interaction || 



! 3 k m 2 - 

V lls — o 2 3 !' 
2 &T 6 ITli 



the interaction of the orbital angular momentum I with the spin s 2 of the central body |7j] 



and the spin-spin interaction 



V L ( 2 ) 



k 

Vss = ^3 [ 3 (sin)(s 2 n) - sis 2 ]. (3) 



Simple symmetry arguments dictate now the form of the spin-orbit interaction for the two 
body problem: 

3 k ?(m 2 _ mi 



V 2ls =2-^l{si+s 2 ). (5) 



As to the spin-spin interaction, it is of the same form ([3D. 
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However, due to the mentioned discrepancy concerning the spin-orbit corrections to the 
gravitational radiation, it turns out expedient to derive explicitly the interactions discussed. 
This is only an elementary generalization of the solution of the corresponding problem for 
the case of a heavy central body, as given in book |J (§106, Problem 4). We will start with 
the two-body Lagrangian including c~ 2 corrections: 

m\v\ m 2 v 2 kmim 2 rriivf + m 2 v 2 

km x m 2 2 2 7/--x /-» -»\/-» ->\i k 2 m x m 2 {m x + m 2 ) 

Here r = r\ — r 2 \ n = f/r; rrii, fi, Vi are the mass, coordinate and velocity, respectively, of 
the zth particle, i = 1,2. 

Let us take the term with v\ in the second line of Eq. We write the velocities of 

individual elements of the top 1 (with mass dm\) in the form 

Vi + UJi x p x , 

where v\ is the velocity of the orbital motion, uii is the angular velocity. The radius-vector 
pi of the element dm\ is counted off the center of mass of the top 1, so that the integral over 
the volume of the top 

Pidmi = 0. (7) 



Due to (0) the first-order term of the expansion in pi/r of the interaction discussed vanishes. 
As to the second-order term, with the obvious definition 



J PimPindrrii — — Si 6 mn 



of the spin si of the top 1, it generates 

3 k m 2 
2 c 2 r 3 mi 



si[r x p x \ 



in the spin-orbit potential. Treating in this way the next term, that with v 2 , in (|(|), we 
completely restore the spin-orbit potential (|^) in the center-of-mass system for the binary, 
where pi = —p 2 = p. The similar procedure applied to the terms with — 7(viV 2 ) — (vifi)(v 2 n) 
in @ leads to the next spin-orbit contribution @, as well as to the spin-spin potential (|3|). 
It should be mentioned that the above expressions for the spin-orbit and spin-spin inter- 



action in the two-body problem were obtained previously in Refs. || [K| from the analysis of 
the scattering amplitude for spin-1/2 particles in the one-graviton-exchange approximation. 
(As to the spin-independent relativistic correction, some terms of this type are missing from 
their expression.) 

An amusing fact is that the obtained spin-orbit and spin-spin interactions are exact ana- 
logues (up to an obvious change of notations) of the corresponding well-known terms in the 
hydrogen atom. We mean the fine and hyperfine structure, the last interaction being induced 
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by the coupling between the nuclear spin and the electron orbital angular momentum and 
spin. (Of course, in our classical approach we cannot reproduce the contact Fermi spin-spin 
interaction with 5(f) ). 

However, the expression for the spin-orbit correction to the acceleration, presented in 
Refs. [|I|, |J, differs from that which can be derived from our formulae Vu s and V^s- The 
discrepancy is due to the difference in the definitions of the center-of-mass coordinate of a 
rotating star. The coordinate X{ advocated and used in Ref. [p]] is related to ours f as follows: 

fi= Xi+ x ^ (8) 

(from now on we put c = 1 in our explicit formulae). The shift by itself is of course a matter 
of convention, but it is in fact our definition which just by construction (see Eq. (|^) and the 
arguments leading to it) corresponds to the true center-of-mass coordinate of a rotating star. 

Still, what is the meaning of the vector x and why is it irrelevant to the problem under 
consideration? The answer can be conveniently formulated with the following example. For 
the free Dirac particle with the Hamiltonian 

H = dp + (3 m 



the operator whose expectation value equals to r, is not r itself, but [11 



x = r + ; — - — — ; E v = \ hp + m z \ S = — a x a . 9 

2E P 2E p (E p + m)p p V ' P 2i [ J v ' 

To lowest nonvanishing orfer in c~ 2 expression (P) reduces to 

1 - - - & 
x = r — — v x s; s = — (10) 

which might prompt indeed substitution @. However, the transition from the exact Dirac 
equation in an external field to its approximate form containing only the first-order correction 
in c~ 2 is performed by means of the Foldy-Wouthuysen (FW) transformation. And under the 
same FW transformation the relativistic operator x (its form for an interacting particle is 
more complicated than (|S|)) goes over into mere f. In other words, in the arising Hamiltonian 
the coordinate of spinning electron has the same meaning f as in the completely nonrelativistic 
case. Nobody makes substitution (§) when treating the spin-orbit interaction in the hydrogen 
atom. 



3. Let us consider now the fully covariant equation of motion for a spinning particle in 
an external field 

D ( DS^ W \ 1 

derived by Papapetrou ||12"| . Here D/Dr means the covariant derivative with respect to the 
proper time; = dx^/dt is the four- velocity; S^ u is the antisymmetric tensor of spin; 
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is the Riemann tensor. We have included as well into this equation the interaction with an 
external electromagnetic field F^ v '. A close analogy between the two terms, electromagnetic 
and gravitational, in the rhs of Eq. flTTl ) was emphasized in Ref . [0 . 

We will use the common definition of the relativistic spin. According to it, the only 
nonvanishing components of the tensor of spin (and the vector of spin) in the particle rest 
frame are the space ones. Transition to an arbitrary frame is performed by a boost. This 
definition guarantees automatically the constraint for spin 



S^Uy = 0. 



(12) 



Due to this constraint, 



DS^ 



n U »=- S " v ^r- ( 13 ) 

Dt Dt 

So, if the electromagnetic field is switched off and terms nonlinear in spin neglected, the 
second term in the lhs of Eq. ( |TT|) should be deleted. Clearly, — ^R^u u S pa is nothing else 
but a covariant expression for the force due to the spin-orbit interaction. In the field created 
by a heavy mass M this term reduces to first order in c~ 2 to 



kM _ _ 
3 — — (V x s 



(nv) n x s — 2n(n[v x s])) 



(14) 



which coincides with the corresponding force from Ref. 
from potential (ffl) is different: 



HT1. However, the force extracted 



kM 
3 — — I v x s 



- (nv) n x s — - n(n[v xs])j. 



(15) 



This discrepancy was pointed out long ago in Ref. [I3j where the force (15|) was derived 
from the scattering amplitude for the Dirac particle. The explanation suggested in Ref. |y| 
for the disagreement is that expression (|I"4]) refers to an extended body and (^5j) to a point 
particle. It does not look satisfactory. For instance, is the proton in a gravitational field 
a point particle or extended body? Obviously, as long as we do not go into details of its 
structure, as long as we do not consider its internal excitations, an extended body can be 
treated clS cL point particle. 

To make the problem even more acute, let us consider another limit, that of vanishing 
gravitational field. Here Eq. ( P"T|) describes a particle with spin, but without magnetic 
moment. Still, its spin interacts with an external electromagnetic field, which to lowest order 
in c -2 should be described by the Thomas interaction 



V t = —s[Exv}. 
2m 



(16) 



This expression can be easily recovered from the well-known results for the spin precession 
(see, e.g., book |14]]) at the vanishing ^-factor, g = 0. When the electric field E is that of a 



point charge — Ze, the force corresponding to interaction (|T6|) is 



Ze 2 



V X s 



— [nv)n X s 
2 v ; 



— n(n[v x s\ 



(17) 
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However, the force obtained in this case from the second term in the lhs of Eq. ( |TTD to lowest 
order in c -2 , is different: 

Ze 2 

— — (v x s — (nv) n x s — 2n(n[v x s\)) . (18) 



The reason of both discrepancies is clear now: f entering expressions (14), (0) is just the 
relativistic coordinate of Eq. (|TT|), it is nothing else but x in the notations of relations (§). 
Therefore, the transition from the fully relativistic Eq. (11) to the c~ 2 approximations to it 



([14]), (|18|), should be accompanied indeed by substitution fllpp . This substitution should be 
performed of course both in the acceleration entering the Newton equation of motion, and in 
the (formally) nonrelativistic force. In this way correct Eqs. QI5]) , ( |P7| ) are restored. 

4. The above consideration of the Papapetrou equation ( |TTD is instructive in one more 
respect. It was pointed out that this equation describes a particle with spin, but without 
magnetic moment. The magnetic moment interaction is well-known to be taken into account 
by the following term in the relativistic Hamiltonian: 



V mm = (19) 



In Ref. it was demonstrated that expression ([H5]) has a close gravitational analogue 

Vqm = ~~ a S P Rnupa (20) 

8m 

which can be called gravitational magnetic moment interaction. In particular, this coupling 
arises in a natural way in wave equations, and the value k = 1 for the constant in it is as 
preferable from the point of view of the high-energy behaviour as g = 2 is in the electromag- 
netic case 0. Both interactions, ([19]) and (^0|), should include in fact some additional terms 
treated in detail in Refs. |15[ (for usual magnetic moment) and [[J. Being certainly of higher 
order in v/c, those terms can be omitted in our treatment of binary stars. 

For the field created by a heavy mass M interaction (^) reduces in lowest, first order in 
c~ 2 to the quadrupole form: 

3kM 

V} = -^rQ S mn n mn n (21) 

where the effective quadrupole moment 

Qmn ( S m^n 77 ^mn^ ) 

m 6 

resembles by its spin dependence the well-known expression from quantum mechanics. For 
the two-body problem under discussion expression ( pl|) generalizes to the following self- 
interaction of spins: 



k ( 1712 mi \ 

V s = k — - I — si m si n H s 2m S2n {3n m n n - 5 mn ). (22) 

2r d \mi mo J 



resembling the usual spin-spin interaction @. 
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Let us compare now the effective quadrupole interaction ( pTj) or (p2| ) with the usual 
quadrupole one. At k ~ 1 interaction (|22|) is of the same order of magnitude as the spin- 
spin one (|3|). Even in the most favourable case when they can become essential, that of two 
extreme Kerr black holes, both interactions are of the c~ 4 order. The star rotation velocity is 
here ~ c, but its radius is close to the gravitational one r g ~ c~ 2 , so that each spin s ~ c~ l 
III]. (The same argument demonstrates that the spin-orbit interaction is of the c~ 3 order [[J.) 
As to the usual quadrupole interaction, it is suppressed by the small value of the quadrupole 
deformation and, according to Ref. [[16], can also manifest itself in the case of two extreme 
Kerr black holes only. 

5. We are going over at last to the spin effects in the gravitational radiation of binary 
stars. In fact, the only essentially new correction to the energy loss obtained by us is that 
due to the spin self-interaction and originating mainly from interaction (1221). Our final result 



for the contribution due to the spin-spin interaction (y) coincides with that presented in 
Refs. j[| Qj. As to the spin-orbit correction, our result for it can be in fact obtained from 
the expression given in [Q, |4]] by merely going back to the simple-minded definition of the 
center-of-mass coordinate advocated by us above. However, in parallel with calculating the 
correction due to the spin self-interaction (|22|), we will present corresponding contributions 
induced by the spin-orbit and spin-spin couplings (|5|), (Q). It serves as an independent 
check of the results presented previously in Refs. JT], || (an arithmetical confirmation only in 
the case of the spin-orbit interaction). 

We start with the well-known expression for the metric perturbation h mn at large distance 
R from the source (see, e.g., ||, §110): 

4/c f 1 R 

ipmn(R, t) = - — J drr mn (r, t- R + fn); ip mn = h mn + - S mn h pp ; n = — . (23) 

The source r mn includes not only the energy-momentum tensor of matter, but generally 
speaking corresponding nonlinearities of the gravitational field itself. It is conserved in the 
sense 

dpTp, = 0. (24) 

As usual, we will be interested in the part of the perturbation if) mn which is orthogonal to 
n and trace- free (otherwise Eq. (^) would look slightly more complicated). It should be 
mentioned here that both expression (p3|) and the c -2 Lagrangian (|6]) are valid under the 
same gauge conditions 



d m hmo 2 d h mm -0, d,h, n 2 <9„V-0- 



One can easily check it by inspecting the corresponding derivations in book ( §§106, 110). 
Neglecting the retardation fn in expression (|23|) we reduce it to the quadrupole formula 

2k f 
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To lowest order inv/c t 00 reduces to rest masses, and the integral gives the usual quadrupole 
moment 



Qr 



3 ' ' ' mi + m 2 
Here new terms in the quadrupole radiation intensity are generated by the spin- dependent 
corrections to the orbit radius r and to the equations of motion used to evaluate time deriva- 
tives of r. In all our discussions of gravitational radiation we restrict to the case of circular 
orbits which is most interesting from the physical point of view [|l|], and much more simple 
as concerns calculations. In this way we get the following relative corrections to the usual 
quadrupole formula: 



1 ls 

ss _ 

9k 



1 



ni\ni2T l 
9 



9 - 



s 



Ii 4:m 1 m 2 r 2 



2niim 2 r 2 
m 2 



3 S\ t S2t — S1S2) 



nil 



(3 s 



+ 



m 2 



(3 s 



2/ 



(26) 
(27) 
(28) 



Here 



?>2k^ni\ni\ni 
5r 5 



is the unperturbed quadrupole intensity, m — nil + m<i\ subscripts Is, ss, s refer to the 
spin-orbit, spin-spin and spin-self-interaction contributions respectively; 



s = si + s 2 ; C 



m 2 
nil 



81 + 



nil 
m 2 



s 2 - 



The expressions for Ij] and Ij 1 have been averaged over the period of rotation. That is why 
they contain the spin components s t orthogonal to the orbit plane. 

The next correction to the qudrupole radiation originates from the terms of the relative 
order c~ 2 in Too- The only spin-dependent contribution here is of the Is type. The same 
procedure which has generated the spin interactions from Lagrangian (P) allows one to extract 
from 



niiv\ 



+ 



m 2 v 2 



2 2 

the following correction to the quadrupole moment 



m 



(29) 



Since this expression will be anyway contracted with the symmetric Q m n, there is no need to 
symmetrize it explicitly. Certainly, correction (|29|) makes a spin-orbit type contribution to 
the radiation intensity only. But we will postpone its calculation for the time being. 

Let us go over now to the retardation effects in radiation. The first-order correction in 
formula ([23]) looks as follows: 



ab 



4k 
~R 



d / dfzr, 



ab ■ 



(30) 
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We have made explicit here that the wave propagates along the z axis, a,b = 1,2. Simple 



transformations based on the continuity equation (24) (see ||, §110) lead to the following 
identity: 

J drr k r mn = -dl J dfr k r m r n T 00 + - d J dfr n (r k T 0m - r m r ok ). (31) 

The first, totally symmetric term in this expression generates the octupole radiation. Being 
spin-independent, it is not of interest to us. 

In the second structure we restrict to the term in r 0m which is of lowest order mv/c, and 
obtain in this way 

dfr n (r k T 0m - r m r ok ) = mi(ri m ri n t;i fc - r lm ri n v lk ) + m 2 (r 2m r 2 nV2k - r 2m r 2n v 2k ). (32) 

With the previous trick we single out in this tensor the spin-dependent terms and arrive at 
the expression which can be presented as 

The second-rank tensor here is 

Jmn 7, ( ^rrJ"n f'rj'm) T ( ^"mCn ^nCm 7T '"C ) (^3) 

2m 4 3 

where 

? _ _£i S2_ 

77l\ 7712 

It is a close analogue of the magnetic quadrupole moment in electrodynamics, one can single 
out in it in an obvious way the contributions of the convection and spin currents. 

The intensity of this gravimagnetic quadrupole radiation can be conveniently calculated 
via the following transformation of the initial structure n k e mn J kmn : 

Tlk&mnJkmn ^Ink^mn^lkJlm ^Vm^lm- (34) 

If we choose the independent components of the polarization tensor as 

en - e 22 e 12 + e 21 



then the dual polarization is 

( e 12 + e 2 \ e n - e 22 , 

&lm — £lnk&mn n k — I n ' o ] " \""-V 



Formally the sum over independent dual polarizations e in ( e; m J; m ) 2 looks exactly the same 
as that over common polarizations e when calculating the usual quadrupole radiation. In the 
present case the intensity equals |1|, [I7| 

jgmq _ ]^ k4 t(3) t(3) /o fi \ 
1 _ 45 J rnn J mn \ ov J 



where the superscript at J denotes the third time derivative. The calculation of these deriva- 
tives is simplified in the present case by the fact that to our accuracy both I and can be 
considered as constant in time. The spin-dependent corrections arising in this way are 

jgmq 

t(s-i); (37) 



I q 12mim,2r 2 

jgmq j 



sits 2 t ~ 7sis 2 ). (38) 



I q 48 mirror 2 

Even a spin-self-interaction correction (somehow missed in Ref. |1]) arises here: 

jgmq ^ 

mi 1 " *' m2 



— ( s u - 7s i) + — ( s 2t - 7s 2 ; 



J 9 96 mirror 2 

Let us consider at last the second-order retardation correction in formula 



(39) 



^ ab = -^\dljdfz'r ab . (40) 

A spin contribution can be produced in it only by velocity-dependent terms in r k riT mn , in 
other words only by the energy-momentum tensor of matter. These terms are of the type 
VrnVn^kfi and generate the following structure 

^ ^"k^m^lnr^r- 

m 

The symmetrizations (k, I) and ( m, n) are implied here. The irreducible part of the third- 
rank tensor r^v m £ r can be omitted at once since it does not interfere in the total intensity with 
the second-rank tensor Q mn . Then, we omit also the structures of the type (rkV m — r m Vk) £ r 
since both orbital angular momentum and spin can be considered constant in time to our 
accuracy. The resulting structure transforms as follows: 

-E / m _ f. _ _ 2 // 

flk^-l^mn^lnr^m \^k<,r ^r<,k) „ ^•k^l^mn^lnr^m^krs^ijs^i^j „ ^■mn^m^nrs' rSs - 

3m 3m 3m 

In other words, this correction to the quadrupole moment is 

$Qmn = ~ \ — VmenrsTris- (41) 

3 m 

Adding up expressions (|29"D and (^TJ) and neglecting again the term (r m v r — r r v m )£ s , we 
obtain the following total spin- dependent correction to the quadrupole moment: 

3Qmn 77 ^m^nrs^rCs- (42) 

3 m 

The corresponding relative correction to the radiation intensity constitutes 

It 2 



I q 3 mirror 2 



It. (43) 
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Now, expressions fl26|), (|37|) and fl4~3"D, taken together, give the following total spin-orbit 
correction: 

4, _ _f(73s + 450 



(44) 



12 mim 2 r 2 

It can be easily checked that the corresponding result of Refs. [|l], |||] would be reconciled with 
this one under the proper definition of the center-of-mass coordinate. 

Adding up expressions (|27j ) and (|38|) , we obtain the result of Refs. |], |J for the spin-spin 
correction: 

649si t s 2t - 223s!S 2 ). (45) 



± ss 



48 mim2r 2 

And at last, the total spin-self-interaction correction, generated by (|28|) and (|39"1), is 



1 



I q 4m 1 m 2 r 2 



24 



• (46) 



As has been mentioned already, at k 
spin-spin one. 



m 2 z4 m x m 2 

1 this new correction is quite comparable to the 
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